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>• ■ Abstract 

00 _ 

OO \ We study pairs of interacting self-avoiding walks {cj^, cj^} on the 3d simple 

cubic lattice. They have a common origin cjq = o^q; ^^^ ^^^ allowed to overlap 

only at the same monomer position along the chain: u\ ^ u"^ for i ^ j, while 

^\ [ u}\ = uf is allowed. The latter overlaps are indeed favored by an energetic gain 

e. 

j5 ■ This is inspired by a model introduced long ago by Poland and Sheraga [J. 

Chem. Phys. 45, 1464 (1966)] for the denaturation transition in DNA where, 

^^ ^ however, self avoidance was not fully taken into account. For both models, 

pH ■ there exists a temperature T^ above which the entropic advantage to open up 

O ■ overcomes the energy gained by forming tightly bound two-stranded structures. 

Numerical simulations of our model indicate that the transition is of first 

order (the energy density is discontinuous), but the analog of the surface tension 

S^ . vanishes and the scaling laws near the transition point are exactly those of a 



a^ 



o 



H \ second order transition with crossover exponent cj) = 1. Numerical and exact 

analytic results show that the transition is second order in modified models 
where the self-avoidance is partially or completely neglected. 



1 Introduction 

The study of the nature of the DNA denaturation is a long standing open problem. 
Experimentally a multistep behaviour in light absorption as a function of the tem- 
perature was observed already in the fifties (see |l[] as a review). This suggested a 
sudden sharp opening of clusters of base pairs in cooperatively melting regions. This 
scenario is reminiscent of the behaviour at a discontinuous first order phase tran- 
sition, in which the system changes its state from a double strand to two molten 
single-stranded chains. Since then, this scenario has been verified and studied in 
great detail [Q. 

Early theoretical attempts to model this transition could not reproduce these 
phenomena. The first attempt with a one-dimensional Ising model in which the two 
states of spin correspond to a open or close state of the base pair, with a favorable cou- 
pling between neighbor pairs that are in the same state 0,0, reproduced a crossover 
between the two different regimes but no thermodynamical transition. 

The first refinement consisted in taking into account the different entropic weights 
of opened bubbles and double stranded segments [Q , since the phase space region that 
two terminally joined (but otherwise free) open strands can explore is bigger than the 
one accessible to a double strand of the same length. 

This model was solved using the entropic weights of self-avoiding loops in refs. 
1^,^]. In this way the self avoidance between bases within the same loop is taken 
into account, but the other mutual excluded volume effects are completely neglected. 
This simplified model displayed a smooth second order transition in two and three 
dimensions. 

These models were of course only very rough caricatures of the true complexity of 
the problem. Even if we believe that microscopic details should be irrelevant for the 
existence and order of the DNA melting transitions, there are a number of aspects 
which one might suspect to be relevant. In addition to self avoidance these include 
the stiffness of DNA, the difference in stiffness between single- and double-stranded 
DNA, the different properties of A-T and C-G pairs, and the helical structure of 
double-stranded DNA. Finally, one should also consider the effect of "wrong" base 
pairings, either between bases of the two different strands or between bases within 
the same strand. 

There seems to exist up to now no model which incorporates all these aspects. 
But there have been recent models where some of them where included, and which 
seem to reproduce the sudden opening of base pairs. The common property in all 
of them is an entropic barrier that favors configurations in which base pairs are far 
apart. 

The 'nonlinear model', introduced in [0,||, assumes that the stacking energy be- 
tween neighboring base pairs depends on whether these pairs are in 'helical' or 'coil' 
states (i.e., whether they are bound in a double string or not). In a helix, this stack- 
ing energy is larger than in a coil. Transfer-integral calculations, molecular dynamics 
simulations and approximate analytical methods [§] pointed out a first order phase 
transition. 

In the same direction goes also a recent model Q in which the helical structure is 
taken seriously. As a result, a mechanical torque which tends to increase or decrease 



the winding becomes a new thermodynaniical variable. A transfer matrix calculation 
p!0| shows that this model exhibits a first order phase transition in the temperature 



- torque plain, analogous to the liquid-gas transition in the temperature - pressure 
plane. 

Finally, according to a recent study [^ , the effect of the heterogeneity in the DNA 



sequence - which amounts to a frozen disorder in the base pair binding strength - 
has no effect on the order of the transition if the model contains no entropic barrier. 
But it gives rise to a multistep energetic landscape if a state dependent stiffness of 
the type considered in [0 is introduced. 

In the present paper we consider a simplified model where all these features are 
disregarded, but - in contrast to the papers mentioned above - excluded volume in- 
teractions are fully incorporated. Our model consists of two interacting self-avoiding 
walks, corresponding to the two single strands, with the same origin on a 3d cubic 
lattice. Each monomer corresponds to a base and is supposed to have its complemen- 
tary at the same position in the other chain. Two monomers with different positions 
in the two chains are not allowed to occupy the same lattice site, whereas the overlap 
of monomers at the same position is favored by an energetic gain e that represents 
the binding energy. Base-pair misalignments are forbidden. We consider the homo- 
geneous case, where all the binding energies are equal. 

In our approach we focus mainly on the two conflicting tendencies of the system: 
the entropic gain due to the larger number of configurations accessible to the two 
open strands system on one hand and the tendency to build energetically favored 
links between the two strands on the other. The necessity to balance these opposite 
tendencies when minimizing the free energy leads to the finite-T phase transition 
between the high temperature swollen phase, and the low temperature phase in which 
finite fractions of the chains overlap. 

2 The model 

Let us define two A^-step chains with the same origin on the 3-d lattice uj^ = 
{ujj, . . . , wJv} and cj^ = {cUq, . . . , uj%} with cjf G Z^ and cJq = c^q = (0, 0, 0). 

The Hamiltonian (or rather Boltzmann weight) which describes a configuration 
(cj^, cj^) of our system is 

^ /-e ^ \ 

CKT = Y[{1 - 5uji,u;^){l - 5^i^^i){l - 5uj2^uj^) exp 77^ ^ (^w^w? • (2-1) 

i^j ' ' ' ' ' ' V^^ i=0 / 

Thermodynamic properties of the system only depend on the reduced variable e = 
—i/KT that we will use in the following. The partition sum can therefore be written 
as 

N 

^^(e) = Ecjv,ne^" (2.2) 

n=0 

where n is the number of contacts, n = caxd{i\ujl = ujf,i > 0}, and c^^n is the 
number of distinct configurations with n contacts (notice that {n)/N is the natural 



order parameter). Alternatively, by introducing a fugacity z we can go over to the 
grand canonical ensemble with partition sum 



cx> oc N 

Giz, e)=J2 ^""Z^ie) = E E c^-n^^' ■ (2-3) 

N=Q N=Q m=0 

While we fix the starting point of the two polymers at a same origin, lJq = lJq, 
we allow the end-points to wander freely in space. This is different from the Poland- 
Sheraga model [|], where also the end-points were forced to coincide, ijj\i = uj%. At 
least in the ideal case, the presence of this constrain does not affect the order of the 
transition and this should be true also when the excluded volume interaction is taken 
into account. However the crossover scaling functions between different regimes are 
not the same in the two models and at the tricritical point, at which the transition 
takes place, different entropic exponents are found. In the excluded volume case there 
is also a topological subtlety: if chains are deformed continuously, non-trivial knots 
are forbidden if the end points never separate. In contrast, in our thermodynamical 
treatment any knots are allowed. But this should not have much influence either. 



3 Approximate treatment 



The system can be represented as a sequence of M superimposed self-avoiding walks 
of length 111, ... , um (^fc > with n = J2k{^k + 1) ~ 1)) which correspond to helical 
domains in DNA where base pairs are bound together, which alternate with M — 1 
bubbles of lengths pi, . . . ,Pm-i {Pk > 1; molten regions). On the lattice, they are 
self- avoiding polygons of length 2pj. The last part consists of two self- avoiding walks 
of lengths 

M M-1 

r({n,p}) = N -J2ni- J2 Pr (3-1) 

All the elements of the sequence must be mutually avoiding, that means that two 
monomers can occupy the same position in the space only if they occupy the same 
position along the chain. 

If this last constraint is neglected, one can factorize the problem and write a 
generating function for the system in terms of the generating functions of single self- 
avoiding walks, polygons, and a pair of self-avoiding walks starting at the same origin. 
This leads to the 'almost unidimensional' phase transitions of 0,|^. 

The partition sum in the flxed-A^ ensemble can be written as 

^Af = E 'S'r({n,p})n^">^» (3-2) 

{n,p} i 



n 



where the sum runs over all possible partitions into helices and bubbles. Wn 
Q<^ri/kT^^^ Y^ _ Q^^ ^^^ Sr{{n,p}) = C2r{{n,p})- Here c„ is the number of self-avoiding 
walks of length n, while C2p is the number of self avoiding polygons of length 2p. This 
partition sum is clearly an upper bound to the true one, since many conflgurations 
are included which are not allowed due to self avoidance. 



The simplified problem can be easily solved in the grand canonical ensemble, i.e. 
by considering the generating function G{z) = J2n=o z^ Zjy. We find 

""^'^ - I - e^ G^^{z)Gv{z) ^'-'^ 

where G\^{z) = Y.n=o^''Wn, Gv{z) = T.n=^z''Vn, Gs{z) = En=oz''Sn. The 
critical behaviour of the system is determined by the singularity of G{z) which is 
closest to the origin, and it can be studied by using the asymptotic forms for the 
number of self-avoiding walks and polygons c^ ~ n^N''~'^, Cn ~ fi^N"''^, where 
fj, is the connectivity constant, which is a lattice dependent quantity. Two cases are 
possible: that the singularity closest to the origin comes from Gs{z) or from vanishing 
of the denominator. Let us focus on the case d = 3. Recent estimates of the critical 
exponents are a = 2 — du = 0.23723(4) (where we used the estimate z/ = 0.58758(7) 
given in 0) and 7 = 1.1575 ± 0.0006 0. 

Let us denote with zw, zy and zs the location of the singularities of the generating 
functions on the real axis. Since 

it appears immediately that z^ = — , zs = ^, while Gv(z) is finite at zy = zs = ^ 

fi 11 fi 

(since one has 2 — a > 1), but it diverges for z > zy- 

In the high temperature regime, i.e. e -^ 0, the singularity of Gs is the first to 
occur and the critical behaviour is that of two self-avoiding walks. This means that 
the system is in the denaturated state, and the corresponding free energy density 
is given by f/KT = log zs = — 21og/x. This is just the entropy density of two self 
avoiding walks with the same origin. 

Since G^yiz) is an increasing function of e, for decreasing temperatures (increasing 
e) the zero Zwv of the denominator in Eq.( |3.3| ) decreases and finally becomes lower 
than Zs- The crossing point corresponds to the melting transition e*. 

It can also be shown that the order of the transition is determinated from the 
singular behaviour of Gv{z) in ■\, namely on the value of the exponent 2 — a [Q. 
Since Gy^^z) is regular in ^ at e = e*, it plays an irrelevant role at the transition point 
and this is independent from the value of 7, i.e. the fact that the helical domanis are 
self-avoiding or ideal does not affect the order of the transition. The free energy for 
e > e* is given by 

/ -21og^ + C(e-e*)i/(i-") + ... for l<2-«<2 
^' \ -21og/i + C(e-e*) + ... for 2-a>2 ^''■^> 



Since 2 — a = 1.76276(6) [0, the approximate solution predicts a second order 
phase transition. 

The main approximation involved in the above treatment is that it neglects ex- 
cluded volume effects which come from the mutual interactions of different bubbles, 
segments and free terminations the one with the others. As already pointed out. 



this means that we overestimate the partition function, and the transition could be 
sharper than predicted by this simple model. 

On the other hand, one can immediately use the above arguments to infer that 
the transition is certainly of second order in the case of interacting random walks, 
since 2 — a = 3/2 there, and there are no excluded volume effects to be taken into 
account. We present in the appendix an exact analytical treatment of the ideal case. 
There we also evaluate numerically the melting value e*, and we study the scaling 
laws whose general structure will be discussed in the next section. 

4 Scaling laws 

Again it is more easy to discuss the problem in the grand canonical ensemble, with 
the fugacity z conjugate to N . The limit iV -^ oo in the monodisperse ensemble 
corresponds to z /" zjye). 

Values z > Zc{e) make only sense after placing the system in a (large but) finite 
volume V. The two polymers are allowed to grow until they fill the volume with a 
finite non-zero density p = 2{N)/V which remains constant in the limit V — > oo. 
In this regime, we actually have two different phases, corresponding respectively to 
molten and undenatured (double-stranded) chains. It is intuitively clear that a non- 
zero density favors the presence of contacts, because contacts in our model reduce the 
volume occupied by the monomer pairs. Therefore one expects that the transition 
point between the molten and double-stranded phases takes place at a lower value of 
the interaction parameter e'^'^'^'^^z) when the fugacity z increases. This phase diagram 
is shown in Fig. |I]. 

Notice that the boundary between the molten dense phase and the short chain 

phase is strictly horizontal, as the attractive interaction plays no role along this 

transition line. Qualitatively, this phase diagram is very similar to that for a polymer 

attached to an adsorbing surface 0,|15| and to surface transitions in magnetic systems 

T6| . But in contrast to the latter, the boundary between the two dense phases is not 



horizontal. 

Using the [z, e) representation, it is clear that the melting transition is a tricritical 
point. Its analogue in magnetic systems with surfaces is the special point |1^. The 



curve z = Zc{e) consists of two parts, a horizontal one for e < e* and a tilted one for 
e > e*. At the melting point {z* , e*), one sees a change of critical behaviour. At finite 
large A^, for e < e* the critical behaviour of a SAW of length 2A^ is observed, while 
for e > e* the system displays a double stranded behaviour. At fixed A^, for e near 
e*, a crossover between the tricritical behaviour and the double stranded one (e > e*) 
or between the tricritical behaviour and the 2A^-SAWs one (e < e*) is observed as e 
tends to e*. The width of the crossover region decreases as A^ increases. 

In the following we shall discuss the scaling laws that would be expected if the 
melting transition is second order. If it is first order, it seems at first not clear whether 
the usual scaling scenario (which is based on the existence of a divergent length scale) 
still holds. We shall see in Sec.|^ that it does hold even then. An analytic study in 
the ideal case (no excluded volume) is shown in the Appendix. 



0.32 
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q = exp(E) 



Figure 1: Phase diagram. On the horizontal axis is plotted the Boltzmann factor 
q = e'^ per bound monomer pair, while the fugacity is plotted vertically. Below the 
continuous line, chains are short with an essentially exponential distribution in chain 
length. At this line, the average chain length diverges. To the left of the triple (and 
tricritical!) point, the line is horizontal (i.e., the critical fugacity is independent of q 
and coincides with the value for normal SAWs). The "molten chains" and "double- 
stranded chains" phases are well defined only for finite volume V , with the chain 
length A^ oc V^. The numerical determination of the phase boundaries is discussed in 
Sec. 6. While the continuous line is very precise (error less than the width of the 
line), the uncertainty of the molten/double-stranded phase boundary is at least as 
big as the symbol size. 



Near a tricritical point, the partition sum is expected to scale as 



G(z,e) 



[z — z 



I -7 



*F((e-e*)/(^* 



(4.1) 



with (f) being called the crossover exponent. The scaling function F{x) is non-singular 
at x = 0, from which follows 



ZN{t*) ~ {i/z*Ym 



(4.2) 



for the scaling exactly at the melting pointQ. 

For e > e* and z — > Zc{e) from below, G{z, e) must scale as the partition sum 



for a (double-stranded) SAW, G ~ const /(zc(e) — 
singularity at some finite value xq where it diverges as 



'. Therefore F{x) must have a 



F{x) 



const 



(xo 



X] 



(4.3) 



^The value of 7* in the ideal case is computed in the appendix. We find 7* — 1 + (see Eq.(7.43)) 



On the other hand if both the extremities are bound together as in the Poland-Sheraga model, it is 
easy to see that the crossover exponent (p does not change, but the absence of Gs{z) in the numerator 
of Eq.(3^) gives a different singular behaviour at e*. It is simple to see, following the same lines as 
in the appendix, that this gives 7* = 4>. 
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and for small e — e* 

z*-z,ie)r^Xo'^^ie-e*y/^;. (4.4) 

Thus the crossover exponent describes how the critical fugacity depends on the 
contact energy in the bound (non-molten) phase. 

Finally, for e < e* and z y z* we must have G ~ const /{z* — z^. This is the case 
if 

F{x) ~ {-xY , X ^ -oo (4.5) 

with some power a, and 

'y* + (j)a = 7. (4.6) 

Performing the Laplace transform one checks easily that eq.( |4.1|) is obtained with 
the ansatz 

Zjv(e) = fi{efN^*'^^ ((e - e*)N^) (4.7) 

with /x(e) = l/zc{e). In order to get the right asymptotics for e 7^ e*, the scaling 
function \l/(x) - which is related to F{x) by a Laplace transform - must scale as 

1 3; I (7-7*)/'/' fQj- X — ^ ±(X). 



The scaling of the energy is obtained by deriving Z^ with respect to e. It is ||I4 



(e-< 


,*\1/<I>- 


-liV 


for e > 


e* 


N't" 






for e = 


e* 


1/(6* 


-e) 




for e < 


e* 



^7v(e) ~ < A^-^ for e = e* . (4.^ 



From this we see that a first order transition is obtained for = 1. The scaling of 
the specific heat is obtained by deriving once more with respect to e. One finds that 
the peak of the specific heat scales as N^'f'~^ and is located at e* + const/N'^. 

One can also look at the system from an extended scaling point of view ]T7|^T^ . We 
define two correlation lengths ^1 and ^2- The second, which we call the geometrical 
correlation length, is identified with the Flory radius of any of the two polymers, 
5 = ((a;]v - a;J)2)V2 = ((^;2^ _ a;2)2^i/2_ j^ follows the scaling law ^2 ~ N" in any 
phase. The first, ^1, is the thermal correlation length. It is defined as the mean 
diameter of the molten 'bubbles'. In the bound phase we expect it to scale with A^ 
and (e — e*) in the same way as the end-to-end distance between the two strands, 
^1 ex Rend = {i<^N ~ ^nYV^"^ ■ ^^ ^^e deuaturatiou transition is second order the 
thermal correlation lenght ^1 inside the bound phase converges in the limit A^ — >■ 00 
to a constant which depends on e, but the value of the constant diverges as e \ e*. 
Exactly at e = e* it should scale as a function of A^ in the same way as ^2 i-e. both 
correlation lengths should be equivalent. 

In contrast, in a usual first order transition we would expect that the thermal 
correlation lenght remains finite as A^ tends to infinity also in the limit e \ e*, but 
we will see that in our system this picture does not hold because of vanishing of a 
surface tension. 

In approaching the transition point from the molten phase -Rend scales as the Flory 
radius and is then not related to ^1. 

Denoting with a subscript T the exponents that govern the scaling laws in the 
thermal parameter e — e*, we can define a thermal correlation length exponent i't by 



assuming ^i ~ (e* — e) '^^ in the limit where we take first N ^ oo and then e ^ e* — 0. 
One has [|f|[l| 

= u/ut . (4.9) 

This can be understood in two ways. First, one can invoke the fact that ^2 ~ ^1 when 
e = e*. Then Eq. (|4.9|) expresses just the fact that d^'^/de and d^'^/dz are related by 
Eq.(El). 

Alternatively, one observes that D = l/u is just the (Hausdorff-) dimension of the 
system, whence the specific heat exponent a^ = 2 — ^ takes the familiar hyperscaling 
form Ot = 2 - Dz/T [0. 

Let us finally discuss histogram methods which have become increasingly popular 
during the last years. They provide expectation values at temperatures different 
from those used in the simulations. In addition, they are used to study finite lattice 



size effects [20, 21]. Near an ordinary temperature-driven critical point, the energy 



distribution scales in a finite spin system of length L as 

Pl{E) ~ L-i/> ({E - {E,))/L"'') (4.10) 



This is different for a first order transition where the distribution has two peaks which 
get increasingly separated when the system size is increased. The minimum between 
the two peaks becomes exponentially deep (with the depth controlled by the surface 
tension between the two phases), and the peaks become arbitrarily sharp in the limit 

L ^ CX3. 

Instead of studying the distribution for fixed lattice size, in the present case it is 
natural to study it for fixed finite N . In view oi E r^ N"^ , one might now expect a 
similar behaviour 

Pn{E) ~ N-^p{E/N'^) . (4.11) 

In Sec. ^ we show that this is indeed true for the melting transition of ordinary 
random walks in dimensions 2 < d < A. There the transition is second order with 
(f) = d/2 — 1, and Eq.(4.11) is correct. But surprisingly Eq.( 4.1l| ) is also correct for 



ordinary random walks in dimensions d > A where = 1 and the transition is first 
order. This can be understood as a vanishing of the analog of the surface tension: the 
cost involved in going over from a molten domain to a bound domain does not increase 
with A^. This is obviously due to the fact that our system is (at least topologically) 
one-dimensional. Indeed we will show in Sec.^ that the same is also true for SAW 
melting in rf = 3. 

5 Simulations 

We use the pruned-enriched Rosenbluth method (PERM) [^, with markovian antic- 
ipation ||2^, which is particularly effective to simulate interacting polymers |2^. In 



the present case the algorithm was implemented in such a way that the two chains 
grow simultaneously (i.e. adding one monomer to the first chain, than to the other. 



again to the first and so onQ). Following the PERM strategy, the whole system grows 



according to the Rosenbluth method |2^ while configurations with very large/very 



small weight are cloned/pruned. The bias used during the Rosenbluth sampling is 

corrected by multiplying the weights of the configurations with the appropriate factor. 

The /c-steps markovian anticipation consists in an additional bias based on the 



statistics of sequences of fc + 1 successive steps [0]. In dimension d, labeling by s = 
0,...2d — l the 2d directions on a hypercubic lattice and by 5 = (s_fc, . . . Sq) = (s, sq) 
a given sequence of k steps ending in sq, one considers the statistical weight PN,m{S) of 
all A'"-step chains which followed the sequence S during the steps N—m—k, . . . , N—m. 
The bias in a fc-step markovian anticipation is then given by 

2d-l 
p(So|s) = PN,mis, So)/ Yl PN,mis, Sq). (5.1) 

This means that a step in the direction sq is chosen more often if the previous expe- 
rience tells that it will be more successful in the far {m steps ahead) future. These 
biases are obviously compensated by a factor oc l/p{so\s), to get a correct sampling. 
In our simulation the weights PN,m{S) were estimated in a preliminary run. Moreover 
we used an ad hoc bias for the present model. When the second chain has to perform 
a growth step and the end of the first one is in a neighboring site, instead of doing a 
blind step multiplying then the weight by a factor e^ if the new contact is formed, we 
favor it choosing the step towards the end of the first chain with probability oc e*^. 

One can further increase the probability of sampling configurations with many 
contacts by favoring growth steps which reduce the end to end distance r = la;^ — ci;||. 
We found that such a bias (which has to go to zero for r -^ oo) can substantially 
enhance the sampling efficiency, but leads occasionally to "glitches" where a disfa- 
vored configuration is encountered nevertheless, with exceptionally large weight. We 
therefore use only the previously described bias. 

As one could expect for a first order transition, fiuctuations near e* are very 
strong, particularly on the probability distribution of the energy. A large part of 
our statistics was collected in order to obtain clear data for P{E) up to large chain 
lengths {N = 3000, i.e. a total of 6000 monomers), and for a very wide range of E. 
This aim was achieved by performing independent runs at several different values of 
the interaction strength and by reweighting results. The errors were evaluated with 
the jackknife method, i.e. by using the fluctuations between independent runs. Since 
we have only few such runs to compare with, the errors should be considered just as 
rough estimates. 

Moreover we made runs up to very large chains lengths (A^ = 8000) in order to 
study the large A^ behaviour of the partition function and of the end to end distance. 
We also performed simulations at flnite density (which will be discussed more in detail 
in the following) in order to locate the molten/doublestranded phase boundary. 

As a last remark, we note that also in the ideal case in 5 dimensions, where it is 



^This can straightforwardly be extended to more than two chains and allows then very efficient 
simulations of star polymers, in particular near to collapse transitions where other methods become 
inefficient. 
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much more simple to get good statistics, the use of a bias is necessary for correctly 
sampling the energy distribution in a reasonable CPU time. 

6 Results and discussion 

6.1 The ideal case 

Before discussing self avoiding walks, we study first the case of ideal random walks, 
but of course weighting configurations with n contacts by a factor exp(e n). The 
study of the ideal system, which is analytically solved in the Appendix, allows us 
both to test the efficiency of the numerical methods and to verify the peculiar first 
order transition predicted in 5 dimensions. We limit the analysis to simple hypercubic 
lattices with (i = 3 and 5. 
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Figure 2: Scaling plots of the average contact number for random walks in Sd, and of 
the corresponding probability distribution at e = e*. Here e* = 1.07726 and = 1/2 
are the exact values (see Appendix). 

6.1.1 d=3 

In 2>d one finds a smooth second order transition when moving along the solid line in 
Fig. |I|, characterized by a value = 1/2 of the crossover exponent. In Fig. ^ we plot 
{n)/N^ as a function of (e — e*)N'^. The behaviour agrees with the expected scaling 
law, i.e. the data are compatible with the exact values of e* and of 0. Furthermore, 
finite size corrections appear to be small. The probability distribution P{n) exactly 
at the critical point, properly rescaled by A^"^ and plotted against n/N'^, is shown 
in Fig. 1^. The perfect data collapse confirms both the validity of the scaling law 
Eq. ( [4.11|) and the efficiency of the numerical method. 

6.1.2 d=5 

Let us now turn on the more intriguing 5d case. It is shown in the Appendix that 
the system undergoes a ffist order transition, since the contacts density {n)/N is 
discontinuous at e* in the thermodynamic limit, but shows scaling with a value = 1 
of the crossover exponent. We present in Fig. ^ data for P{n) at the analytically 
calculated critical value e*. This plot is completely analogous to Fig. ^b, but uses 
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the exact value = 1 (see appendix). We see now definitely larger corrections to 
scaling. Also, the curves are slightly cap-convex which shows that {n)/N is not strictly 
discontinuous at e = e*. But this is obviously a finite-size effect. For N = 1000, the 
maximum of P{n, e) jumps from n/N = to n/N = c* ~ 0.4 when e is increased 
from e* to l.OOOSe. In the limit A^ — > oo, on the basis of our numerical results we can 
conjecture that P{n, e*) is fiat between n/N = and n/N = c*. 

On the other hand, it is obvious that P{n) does not have the double-peak structure 
familiar from usual first order transitions. In these usual cases, the valley between 
the peaks is due to surface tension: States with the order parameter between the two 
peak values contain more than one domain and are therefore suppressed by surface 
tension. In the present case, there is no penalty for a transition between a molten and 
a double-stranded domain, explaining why all configurations with < n/N < c* can 
be equally much populated (for magnetic polymers, a similar scenario was proposed 
recently in [^). Notice that scaling as formulated in the previous section would 
not hold for usual first order transitions, since the valley between the peaks becomes 
exponentially deep in the thermodynamic limit, P ~ exp(— crL'^~^) for systems of 
linear size L. Formally, this can be reconciled with the present case by noticing that 
our polymers have topological dimension equal to one. 

Related to this are strong fiuctuations in the separation of the two chains in the 
double-stranded phase. Indeed, as proven in the appendix, the average the end-to- 
end diverges as ^i ~ (e — e*)~^/^ as the critical point is approached from e > e*. 
This shows that the 'thermal' correlation length exponent is z/y = 1/2 = 1/D and 
ax = 2- Dut = 2-0-1 = 1. 
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Figure 3: Probability distribution of contact numbers for random walks in bd at the 
exact critical value e* = 2.00115. 

In Fig.^ we plot the contact number density (left panel) and the specific heat 
(right panel). 



C = AT- II = AT-e^^ = N-h-[{n-) - (n)^]. 



(6.1) 



against e. In Fig.^, curves for different chain lengths considered cross nearly at the 
same e value which is also the value the specific heat maximum is moving towards. 



12 



O.bb 
0.5 

0.45 


N=250 ^ 
- N=500 — 

N.750 

-N=1000 


' — ' 


0.4 






0.35 






0.3 






0.25 






0.2 






0.15 






0.1 


-..•*•*'** 




0.05 

n 


liianssss 


•''•'° 



N=250 > 
N.500 ■ 
N.750 ■ 
N.I 000 




2 2.01 2.02 2.03 



Figure 4: Average contact number and specific heat for random walks in 5d. 



Witliin tlie expected uncertainty, it agrees witli tlie predicted e*. But wlien trying 
to make tlie data collapse by plotting them against (e — e*)N'l', we see again strong 
corrections to scaling. Indeed, such fits using [n) /N"^ would give (f) ~ 0.93 instead of 
= 1. Another estimate for the crossover exponent can be recovered from the way 
in which the maximum of the specific heat moves toward e* as the length of the walk 
increases. Since in the crossover region the specific heat follows the crossover scaling 
the position of the maximum for N finite is shifted by a term proportional to A^"*^. 
> 1 by a small amount. These deviations from perfect scaling give us a hint of 
what we have to expect when going now over to self avoiding walks. 

6.2 Self Avoiding Walks, d = 3 

6.2.1 Scaling of P{n) and Properties Derived from it 

In Fig. 1^ we show log P(n,e*), where e* = 1.3413 ± 0.0004 is our best estimate of 
the critical e value. This distribution is clearly not convex, in contrast to the case of 
5-dimensional ideal random walks studied in the previous subsection. Instead, there 
is a peak at n = 0. Due to it, the maximum of P(n, e) jumps discontinuously when e 
passes through e* (see Fig. |^). Apart from this, the situation is very similar. We see 
again substantial corrections to scaling, but it seems quite clear that scaling works 
with = 1. In particular, the depth of the valley between the peak at n = and the 
shoulder at n/N = c* ~ 0.5 does not increase with A^. And, what is more important, 
the value of c* does not substantially decrease with A^. This is our strongest numerical 
evidence for the transition to be of first order. It should be noted that here we are 
considering up to quite long chains (A^ = 3000) and that the data for A^ = 2500 and 
for A^ = 3000 are nearly indistinguishable. 

Deviations from the scaling behavior are seen mostly for large n/N . There, the 
distribution becomes increasingly steeper with A^. This was seen also for ideal random 
walks in (i = 5. It is indeed easy to understand. At large n, we expect P{n) ~ 
e^" (N/n)"' (1 — n/N)"'~ which does not follow our scaling law but is in qualitative 
agreement with our data. 

Accordingly, also the scaling of {n)/N becomes poor for large e, as seen from 
Fig. 1^. If we would try to optimize this scaling plot, we would find ~ 0.94. But 
we know that this would be wrong since this would put too much emphasis on the 
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Figure 5: Probability distribution of the contact numbers for SAWs in 3d at the 
estimated critical value e* = 1.3413(4). 
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Figure 6: Probability distribution of the contact numbers for SAWs in 3d at e 
.999e* and e = l.OOle*. 



region n/N ^ 1 in the contact number distribution which we know to be not scaling. 

We face a similar situation when looking at the specific heat. As seen in Fig. |^, the 
height of the maximum (see Fig. ^) increases roughly oc N, which would correspond 
to = 1, but a least squares fit of our data shows an effective exponent < 1. Since 
we are working with N < 3000, in a regime which is far from the critical one, we 
expect that our estimate is biased by the presence of corrections to scaling. Anyway 
no clear trend in the computed exponent is obtained from fits of data with cuts 
corresponding to increasing values of N (see table (Q)). All the obtained values are 
compatible in the statistical errors, but a x^ analysis indicates that our errors could 
be overestimated. Nevertheless the central value of the effective exponent increases 
with the increasing of Nmin- This could indicate the presence of corrections to scaling. 
The last value obtained for A^ = 1500 is compatible with 0=1. 

On the other hand, fitting the shift of the position of the maximum, which is 
shown in Fig. ^b, one would get = 1.31(14). In this case the large statistical errors 
on our data do not allow any guess on the correction to scahng effects. 

In summary, these results show that the melting transition for interacting SAWs in 
3 dimensions is first order, while the analogous transition for interacting ideal walks in 
3 dimensions is second order. This agrees with our expectation that excluded volume 
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Figure 7: The average contact number plotted against e (left panel) and against 
(e — e*)A^ (right panel). 
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Figure 8: The specific heat. 



effects should make the transition sharper. We should add that we also performed 
simulations of a third version of the model, in which the two polymers were self- but 
not mutually avoiding (data not shown). While two monomers of the same chain 
were not allowed to occupy the same site, two monomers from different chains were 
allowed to do so, and this contributed to the energy if and only if the indices of the 
monomers were the same. For this intermediate model we found a P{n) very similar 
to Fig. ^, but with a less pronounced peak at n = and a very small non-convex 
region. Therefore we are not sure whether the transition in this model is first or 
second order. 

6.2.2 The Thermal Correlation Length 

Up to now we only looked at the scaling behavior of P{n) and at quantities which are 
related to it straightforwardly. We mentioned already in Sec. § several scaling laws 
which are less directly related. 

The most interesting concern the thermal correlation length (see Sec. 4). In the 
double-stranded phase (e > e*) we can identify it with the rms. distance Rend between 
the endpoints of the two chains. Numerical results for the latter (both for e > e* and 
for e < e*) are shown in Fig. ^. For e > e*, -Rend tends for A^ -^ oo to a constant 
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A^min 


0cff ± A0cff 


x' 


DF 


500 
1000 
1500 


0.920 ±0.019 

0.925 ±0.055 

0.98 ±0.15 


0.042 
0.039 
0.0015 


3 
2 

1 



Table 1: Effective crossover exponent 0eflf from the fits of tfie maximum of tlie specific 
fieat as a function of minimum value of Nmm considered in the fit. 
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Figure 9: The specific heat maximum and the corresponding e value for SAW's in 3(i. 



which diverges when the tricritical point is approached, showing that there is indeed a 
divergent thermal correlation length which is independent of the system size for large 
systems. Exactly at the tricritical point we find i?end ~ N" , showing that there the 
thermal and the "geometrical" correlation length (the Flory radius) coincide. The 
latter is also known from polymer adsorption to a wall, e.g. |]T^[T^. The divergence 
of ^1 = limTv^oo -Rend for c — *> e* + is shown in Fig. [Tl|. We see that the thermal 
correlation length exponent Vt-, defined by ^i ~ (e — e*)^'^^, agrees with the Flory 
exponent v (the geometrical correlation length exponent). Since = v/vf [J18|,|19|, 
we find again = 1. 

For e < e* we also have Rend ~ N" . But in this regime, i?end does not scale as 



the thermal correlation length. Instead, we can identify ^i with the average diameter 
of small molten 'bubbles' (which we did not measure since our algorithm would give 
very large errors). Since SAWs in 3 dimensions are not recurrent, large bubbles do 
not occur in the molten phase, and ^i is finite and decreases with decreasing e. 
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Figure 10: Average squared end-to-end distance Rj^^ for various values of g = e% 
plotted against A^ on a double logarithmic scale. Since all curves are based on inde- 
pendent runs, their typical fluctuations relative to each other indicate the order of 
magnitude of their statistical errors. 



Finally, end-to-end distance distributions are shown in Fig. [T^. For e < e* they co- 
incide for large A^ with the end-to-end distance distributions of non-interacting SAWs 
of length 2N, except for a region of very small distances which becomes irrelevant in 
the limit N ^ oo. More precisely, if we denote with cn{R) the number of A^-step 
walks whose end-points are at a distance R = \R\ apart, we have 



fwix) = -Rend-Pw(-R) - -^endZ; - 



Y^RCNiR) 



/W 



-ON 



T-A 



(6.2) 



where x = {2dy^^R/Rend, f{x) is a universal function, and A is a correction to 
scaling exponent. As shown in [^, f{x) is well approximated by a phenomeno logical 
representation given first by Mc Kenzie and Moore [^,|2^. A comparison with the 
latter is shown in Fig. |r^. For small x the attraction between monomers is felt and 
fwix) is larger than for ordinary SAWs, but this effect disappears for A^ ^ cxd, as 
long as e is strictly smaller than e*. The transition between the regimes e > e* and 
e < e* is not through a double-peaked distribution as one might expect for a first 
order transition, but there is an approximately flat plateau at e ~ e*. 



6.2.3 The molten / double stranded phase transition 

The boundary between the two dense phases in Fig. 1 (crosses) was obtained by taking 
finite lattices of size L^ with periodic boundary conditions. During the simulation we 
measured only the energy (number n of contacts) and the partition sum. The phase 
transition is seen as a rapid increase of (n) when the monomer density p = N/L^ 
increases above a value pi (e) . Let us define Zi (e) such that 



dN 
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Ar=pi{€)L3 
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(6.3) 



17 



100 



T 10 



••^, const / ( E - E* f'^^ 



0.001 



0.01 



0.1 



Figure 11: Thermal correlation length ^i = limjv_^oo -Rend plotted against e — e* on a 
log- log plot. For large values of e — e*, statistical errors are much smaller than the 
symbols. For e — e* < 0.01 we cannot give exact errors, but rough estimates can be 



obtained by comparing with Fig. |T^. The dashed line has slope -0.59, corresponding 

to Ut = v- 



If the transition were second order, zi{e) would be the critical fugacity at the consid- 
ered value of e. But our simulations show clearly that the transition is first order. In 
this case, there is a second value of the density /52(e) > Pi(e) such that the system is 
in the double-stranded phase when p > p2, in the molten phase when p < p2, and in 
a mixed phase in between. Equation ( |6.3|) holds in the entire interval pi < p < p2- 

For most values of e we used L = 32, only for e > 1.3 {q > 3.7) we used larger 
lattices of size up to 64^. Using the algorithm as described in Sec. 4 we were able to 
see the first threshold pi, but not the second one. Obviously, our algorithm is not 
efficient enough. Also, the increase of (n) was too slow for p > pi. On the other hand, 
lattices of different sizes gave roughly consistent values of pi. The algorithm became 
much more efficient when we based the "population control" (copying and pruning) 
on weight factors calculated with a larger value of e than that used for evaluating 
average values (we used mostly e = e*, but slightly larger values also worked well). 
This is easily understood. When the density is high enough that double-stranded 
configurations would be favored over molten ones, our ensemble just has no such 
configurations, with overwhelming probability. Thus the parts of the chains grown 
at ambient densities p > pi are correct, but the older parts are wrong and have 
virtually no chance to be corrected by re-growing. This is no longer so when we 
base the population control on a higher value of e. Then chains with larger n are 
favored already from the very beginning. As long as p < pi, they will not contribute 
significantly because they have a tiny weight. But at p > pi, their typical weight 
will become larger than that of the completely molten chains, and they will start 
to contribute. We found that (n) indeed increased very rapidly in a very narrow 
density interval of a few percent, when using this improved algorithm. This still did 
not allow us to measure reliably the difference p2 — Pi, but it convinced us of the 
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Figure 12: Distributions of end-to-end distances at A^ 
q = e''. Normalization is arbitrary. 
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Figure 13: Comparison of the scaling functions /Ar(a;) for A^ = 700, e = 1.19 and 
N = 3000, e = 1.32 with the phenomenological prediction for f{x) of [^,^ . 



correctness of the scenario (in particular of the first order of the transition), and it 
gave uncertainties roughly as large as the symbol sizes in Fig. 1. 



6.2.4 The short chains / double stranded phase transition 

In contrast to the above, the determination of the phase boundary between the phase 
having only short chains and the dense phase for e > e* was straightforward. Values 
Zc{e) were obtained by plotting log Ziy{e) — (7 — 1) log N + N log z against log N and 
changing z until this became horizontal for large A^. Several such curves are drawn in 



Fig. 0. Notice that we cannot expect them to be flat at small values of A^, since the 
behavior Z^ ~ fi'^N'^~^ is expected only when the length of the (double-stranded) 
chain is much larger than the typical size of a molten 'bubble' and this makes difficult 
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Figure 14: Logarithms of partition sums at different values of g = e^, after subtracting 
suitable terms so that they become flat for N -^ oo and e > e*. The straight dashed 
line has slope 7* — 7 = 0.93. 




Figure 15: Phase boundary between the short chain and double-stranded phases. 



the determination of the asymptotic regime. But this does not affect the uncertainty 
of Zc{e) = l//i in a dramatic way, since Zj^j is extremely sensitive to even tiny changes 
of /i. As a result we cannot give formal error estimates, but they are definitely smaller 
than ±0.00005. This number refers to e ~ e* where errors are largest, for e ^ e* our 
estimates are indeed much more precise. A blow up of this phase boundary is shown 
in Fig. [1^. From that we verify that the boundary terminates at {e*,z*) with finite 
non-zero slope, showing again that 0=1. For very large values of e we see also that 
[/i(e)]^e'^ —>■ /i(0) as we should expect for a very tightly bound double-stranded chain. 
From the curves with e ~ e* in Fig. |14| we can read off the exponent 7* controlling 
the partition sum exactly at the tricritical point. In view of the substantial finite size 
corrections at small A^ and the uncertainties at large A^ due to statistical errors and 
the uncertainty in the exact value of e*, we obtain a rather crude estimate 
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2.09 ±0.1 



(6.4) 
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Unfortunately, we do not have any prediction for 7* to compare this with. But for 
ideal random walks we have a^w = ~1 (see Appendix), where a was defined in Sec. 4 
and is related to 7* by Eq.([4.6|). Trying a = Orw as a first guess, we would predict 
7* = 2.16 in surprisingly good agreement with our direct estimate. 



Finally, we can read off from Fig. ^ values of the function A{e) in the scaling 
ansatz 

Z^(e) ^ A{e)i,{ef N^~^ . (6.5) 

According to Eq.(|4.7|), this function should scale as A{t) ~ (e — e*)'''"'^*. This was 
reasonably well satisfied. 



7 Conclusion 

We studied a simplified model for DNA denaturation. It can be considered as a 
lattice realization of the Poland- Sheraga model @, 0, however taking into account 
correctly all excluded volume effects. The melting transition in this model is simply 
related to the balance between the entropic gain of the two DNA strands when being 
independent, to be compared with the energy gain in configurations where they are 
tightly bound together. 

The numerical results show that excluded volume effects are relevant in this tran- 
sition. They change the transition from being second order (without excluded volume 
effects) to being first order. In spite of the latter, we find a divergent length scale, 
and scaling relations typical for a (tri-) critical point. We explain this by the absence 
of a significant surface tension between the molten and bound phases. 

The present study seems the first work where excluded volume effects are taken 
into account systematically in DNA melting. But we have neglected a number of 
other features which might be equally or even more important. We plan to take some 
of these into account in future work. This includes in particular heterogeneity of the 
chain due to the quenched sequence of DNA bases. 

Acknowledgements: We thank Sergio Caracciolo, Walter Nadler, Andrea Pelis- 
setto and Lothar Schafer for discussions and for carefully reading the manuscript. 

Appendix: Exact results for interacting ideal ran- 
dom walks 

In this appendix we provide an analytic solution for our model in the ideal case, 
i.e. when the excluded volume interaction is fully neglected, on (hyper-)cubic d- 
dimensional lattices. In this case the generating function can be easily computed and 
the order of the transition can be easily determinated in any dimension. Let us define 
C2JV (x) as the unnormalized weight of the configuration of two A^-steps random walks 
{a;^,ci;^} starting both at the origin and ending at a distance vector x = (ji\ — ui\ 
apart. A step towards the same lattice position (x = 0) is favoured by an energetic 
gain — e, corresponding to a Boltzmann weight e'^ (the energy is here expressed in kT 
units). We can then write the following recursion relation 
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+ c^jv-2(x-e,-e,)][l + (e^-l)(5,,o] • (7.1) 

Equation ( |7. 1| ) suggests another interpretation. We can choose the origin of our 
system to be sohdal at any time step with the position of one of the two walkers. 
In this way the problem is mapped in the one of a single random walker moving at 
"double speed" (i.e. making two random steps on the lattice at each time step). We 
use here the obvious fact that a walk on a hypercubic lattice can return to the origin 
only after an even number of steps (or that the sublattice with even parity is still a 
hypercubic lattice). This would not be true, e.g., for the triangular lattice. 

If we denote with c^(x) the weight of the single walker configuration, the recursion 
relation ( [7.1[) is equivalent to 

c^(x) = E[c^-i(x + e.) + c^_i(x - e,)][l + (e^ - l)5.,o] , (7.2) 

but because of the rescaling of time only even number of steps correspond to the 
original system. 

Introducing the generating function in the grand canonical ensemble by a Laplace 
transform 

oo 

G^(x,^)=5:c^(x)^^, (7.3) 

N=0 

one can write the Fourier transform G{c[, z) in terms of the free propagator of the 
Gaussian model on the lattice 

^(q'^) = ^^' (7-4) 

where m^ = ^^^^ and Qi = 2 sin ^ , in the following way 

G^(q, .) = ^^^^ , (7.5) 

where the integration is done on the first Brillouin zone. The knowledge of the singular 
behaviour in z in the small-g region in Eq. (|7.5| ), makes it possible to determinate the 
critical properties of the system in the monodisperse ensemble by an inverse Laplace 
transform. 

For instance the partition sum Zj^ = X]xC2Af(x), the mean value of powers of 
components of the square distance from the origin (or in the original system of the 

square distance of the two walkers) (x^™") = "" ^^ ^^ — , and the average number 

of contacts {n) = — ^^ ^^ can be obtained by inverse transforms from the following 
quantities: 

oo 

mo,z) = EE^''c^(x) (7-6) 
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i)'"^^!S^k=o = EE-'^^^c^W 



dqi 



(7.7) 



X Ar=0 



(7.9) 

Because of the isotropy of the system, we write (x^*") = (x^"*). If there were no 
interaction (e = 0), the critical behaviour of the system would arise in the limit 
ml — >■ 0, corresponding to z* = l/2d. We are interested in changes with respect to 
this critical behaviour due to the interaction with the origin. 

It is immediately clear that a critical behaviour which is different from that of free 
random walks case can appear in two cases: If the denominator in eq. ( [7.5|) either 
vanishes or diverges at a z < z*. 

The properties of the denominator are deeply connected to the ones of the integral 
of the free propagator, H{z). 

It is easy to see that 

f {miy^^~^ for d<2, 
H{z) = j [dq]''D{ci,z) ~ <^ logm^ for rf = 2, (7.10) 

[ finite, for d > 2 

so that the case in which the critical behaviour is modified because of divergencies in 
the denominator at z = z* can arise only for d < 2. 

A critical behaviour which characterize a collapsed (double stranded) phase arises 
when the relation 

can be satisfied for Zc < z*, while an unstable fixed point is reached when Eq. ( [7.11[ ) is 



satisfied for Zc = z* and the singularity due to the interaction with the origin merges 
with the one of the free propagator. 

We observe that H{z)/z is an increasing function of z and that lim^^o H{z)/z = 
1, Qfor this reason the above equation has no solution for e < 0, the case of repulsive 
interaction. But in the region e > one can find a critical value e* for which Eq. 
( [7. lip is satisfied for Zc(e*) = z*, with z* defined above. For e > e* the position of 
Zc{e) in Eq. ( |7.11| ) moves closer to the origin, and lime^oo Zc{e) = 0. The point (e*, z*) 
is tricritical in the sense of having co-dimension 2 (two control parameters have to 
be adjusted to obtain it). Eq. ( |7.11| ) gives the critical value e* on the cubic lattice 
which was used in Sec.flO 



.• = -log(l-2<ii/(l)), (7.12) 

Three kinds of critical behaviour in z can be identified in any dimension d, the 
one for (e < e*,z = z*) in which the denominator in (|7.5|) does not vanish (molten 
phase), the one in the double stranded phase governed by (e > e*, 2; = Zc{e)), and the 

■^ The assertion follows immediately from H{z)/z = J2'n=o -2^c|^''(0) 
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one at the tricritical point (e*, z*), but the behaviour in the different regimes depends 
on the diniensionahty. 

• d<2 

From geometrical arguments it is easy to see that when the dimension of the 
space is smaller than the Hausdorff dimension of the walk {dn = ^jv = 2) 
the probability that the walker intersects a given point is finite. In this case 
the interaction with the origin results in a relevant perturbation that brings 
the system out of the universality class of the unperturbed case, as it appears 
from the critical behaviour for ttiq ^ of H{z). The case d = 2 in which the 
dimension of the space is equal to the dimension of the walk corresponds to a 
marginal interaction and logarithmic corrections appear in the scaling laws. 



Eq. ( |7.11| ) has a solution Zc{e) < z* for all values of e > 0, and Zc{e) — » 



when e — >■ +oo. This means that in the attractive regime the free propagator 
stays always finite when the denominator vanishes. In other words there is 
no transition between two different regimes when the temperature is changed if 
the interaction is attractive. The system is in the double stranded phase (or 
collapsed onto the origin in the other view). On the other hand we can identify 
a phase transition when the sign of the interaction changes, namely at e* = 0, 
and we will show that it is a second order one. 

d>2 

In dimension d > 2 the presence of a zero-dimensional region that couples with 
the random walk introduces an irrelevant operator. The integral H{z) takes a 
finite value and its expansion around the singularity z* of the free propagator 



has the form pO | 



H{z) = J2 Anml^ + C (ml) m^^^ m^ = (1 - z/z*)/z . (7.13) 

n=0 

where An are suitable constants and C {m"^) is a function of tuq which is finite 
for ?tiq -^ for d odd and diverges logarithmically for d even. 

The leading term in H{z) — H{z*) for dimension rf < 4 is given by the non- 
analytic one C (ml) mg"^, it is m^ log(mQ) for d = A, and for d > 4 it is given 
by the analytic term tUq. We will show that this gives rise to a change in the 
order of the transition when d = 4 is crossed. 

Molten phase: 

In d > 2, for e < e* the critical behaviour in the generating function is governed by 
the free propagator -D(q, z) whose behaviour around z* gives rise to the usual critical 
exponents of the random walk, but with amplitudes depending on e. One can rewrite 
near to the critical point 

z*e^ + (1 — e^jH{z*) 
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and, from the small q critical behaviour of G^(q, z) and of its derivatives with respect 
to q and e, one can obtain the following scaling forms for large fixed A^: 
partition sum: 

free energy per base pair: 

/= lim ^\ogZl^ = \ogz*^ (7.16) 

W— too iV 

base pair separation moment: 

(x^"^) ^ (2m - l)!!2*"4'^Ar™ (7.17) 

average number of contacts: 

In the limit e ^ e*~ the mean number of contacts (n) and the amplitude of Z\^ 
diverges 

^2jv ~ ,,,,/\ r(e*-e)-iz*-'^ (7.19) 

(n) ^ (e*-e)-i (7.20) 

(7.21) 

For (i < 2, instead, the singular behaviour of Hi^z^ near z* plays a role in the molten 
phase. For e < e* = 0, near to the singular region one has 

G-(q,.)^ ^^^'"| (7.22) 

and the above scaling laws become 

^e /(l-e^)-'z*-2^iV^/2-i forrf=l .„_. 

^2^ ^ \{\-er\\ogNY^z*-^'' forrf = 2 ^^'^^^ 

/ = log(^*') (7.24) 

(n) ~ p^ (7.25) 

(a;2™) ~ A^™ (7.26) 

(7.27) 

from which the behaviour for e ^ 0^ can be easily recovered. 



25 



Double-stranded phase: 

Let us study now the collapsed phase e > e*. The denominator determines the critical 
behaviour, since it vanishes for a Zc{e) < z*, where the small q behaviour of D{q^,z) 
is not singular. Near to the critical region the generating function can be rewritten 
as 

K{e)Zc{e) (1 - z/zc{e)) 

where K{e) = (e" — l)^\zc{e) — e^ is positive. The derivative with respect to e of ([7l5| ) 
in the same limit reads 

dmci,z) ^^ e^{H{z,{e))-z,{e)) 

de z,{tYK\t){l~z/z,{t)f ^^' ' ^ ' ' 

In the same way as the previous relations for the observables have been recovered, 
one finds that 

"^27v ~ T^i \ /I / \ / *\ (7.30) 

K[t) (1 - z^{t)/z*) 

f = log(z,(e)2) (7.31) 

{x'n ^ {2m-l)\\m\- y., .^ (7.32) 

(1 - z^{e)/z*) 

It is interesting to study the limit e -^ e*^ , whose comparison with e ^ e*~ gives 
the order of the transition. The behaviour of the denominator near to e* and the way 
in which Zc{e) approaches e* will be crucial. 

We have to distinguish different cases with respect to the dimensionality of the 
space. 

Let us start with d > 2. From Eq. ( |7.1(J| ) it follows that 

f (1 - z/z*)'^^^'^ for 2 < rf < 4, 

ze'* + (l-e'')H{z)r^(l-e'')l(l-z/z*)\og{l-z/z*) for rf = 4, (7.34) 

[ (1 - z/z*) for rf > 4 

From this relation it is a simple matter to get the behaviour of z* — Zc{e) and K{e) 
for e ~ e* 



2 



z - zJe ~ 



(e-e*)^^ for2<d<4, 

(e-e*)(log(e-e*))-' for rf = 4, (7.35) 

(e - e*) for d > 4. 

(e-e*)^ for2<rf<4, 
^(e) ~ { log (e - e*) for d = A, (7-36) 

constant for d > A. 

The shift exponent ip, defined by 

[e) = z* (l - k {e - e*)'^^) , (7.37) 
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Zr 



with k constant, can be read from the above relations ( [7.35| ) as a function of d. 
Introducing the crossover exponent 



^ for 2 < rf < 4, 
1 for d > 4 



(7.38) 



the above relations can be rewritten in a more compact form, namely z* — Zc{e) ~ 

(e — e*) ' , K{e) ~ (e — e*)^, both with the appropriate logarithmic correction in 
d = A. One can observe that (p coincides with the shift exponent ip. It follows that 
for e ~ e*, 



zAe 



-2N 



^2N 
f 

(n) 



log (z*2) + a 

(2m-l)!!m! 
'H(z*) 



(3* 



{e-e*f'^ forrf^4 

(e-e*)[log(e-e*)]"^ for rf = 4 ' 

(e - e*)-"/'^ for rf ^ 4 

[log(e-e*)]'"(e-e*)-" for rf = 4 

jyf(e-e*)^ forrf^4 

l[log(e-e*)]"^ forrf = 4 



(7.39) 
(7.40) 

(7.41) 

(7.42) 



From (|7.13| ) one sees immediately that the free energy has a discontinuity in its first 
derivative if o? > 4, i.e. the transition is a first order one. 

On the other hand a different critical behaviour is expected at the critical point 
e = e*. In this special point the denominator vanishes exactly at Zc{e) = z*, the 
point in which the free propagator has the small q singularity. The behaviour of the 
denominator of G"* {q, z) , its derivatives with respect to q and d^G^{(\, z)\e* for z near 
to z* can be read from the asymptotic expansion in (|7.34| ). 

Using the same exponent (f) defined above, the following scaling laws hold 



'2N 



{n) 



^^-2N^4, for rf ^ 4 

z*-^^N'*'[\ogN]'^ forrf = 4 



(2m- l)!!m!r(l + 0)r- 

r(m + 1 + 0) 
{N't> for rf ^ 4 



(7.43) 
(7.44) 
(7.45) 



The case d <2 can be handled in similar way knowing the behaviour of H{z) near 
to z*. The way in which Zc{e) approaches z* and K{e) diverges when e — > is given 
by 



^c(e)) 
K{e) 



e2-d ioT d = 1 

g-a/6 for ci = 2 

2 

e 2-d ioT d = 1 

e^/'e for rf = 2 



(7.46) 
(7.47) 
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with a constant. It follows that the above observables in the limit e — > 0+ behave as 



Taking the limit e — > 0+ in Eq. ( |7.49| ) and e ^ 0^ in Eq. ( |7.24| ) it is a simple matter 
to see that the transition is a second order one. 

It's a simple matter to see what is the behaviour at the tricritical point. G'^(q, z) 
becomes the free generating function and the scaling laws that appear are the usual 
ones. The number of contacts can be simply recovered from the derivative of (|7.5| ) 
with respect to e setting e = 0. This gives the well known result 

^"^^1 log AT ford = 2 ^'-''^ 
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